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Abstract 

qq We summarise some basic issues relevant to the optimisation and design of space-based experiments for 

the observation of the Extensive Air Showers produced by Ultra-High Energy Cosmic Particles interacting 
with the atmosphere. A number of basic relations is derived and discussed with a twofold goal: defining 

O^l requirements for the experimental apparatus and estimating the exptected performance. 
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1 Introduction 

Ultra-High Energy Cosmic Particles (UHECP) with energies in excess of « 10 19 eV hit the Earth with a very 
low flux of about one particle • km" 2 • sr _1 ■ millennium - 1 , for particles with energy E > 10 20 eV pQ. 

The observation of UHECP and the interpretation of the related phenomenology is one of the most challenging 
topics of modern High-Energy Astro-Particle Physics. Direct detection is impossible at these energies, due to the 
exceedingly low flux, but UHECP can be detected by observing the Extensive Air Showers (EAS) produced by 
the interaction of the primary particle with the Earth atmosphere. See [2|3j|4], and references therein, for recent 
reviews on these topics. 

The ground-based Pierre Auger Observatory (PAO) [5] is currently taking data: its south site in Argentina 
and its forthcoming north site in the US will provide in the next few years, a clear understanding of many 
important topics [5J [T]. However it is likely that the next generation of experiments for the study of UHECP, 
after PAO, will be space-based, in order to increase the event statistics by exploiting the larger instantaneous 
geometrical aperture which can be obtained by a space observatory with respect to ground based experiments. 

The aim of this paper is to discuss and summarise a few key issues relevant to the design and optimisation of 
space-based experiments for the observation of UHECP. Analytical and semi-analytical relations will be presented 
and discussed, in order to roughly set the scenario for space-based observation of EAS. 

A first discussion of these arguments can be found in 7 , which was a starting point for many of the results 
derived here. 

The results we present are a basic input to the design and optimisation of the experimental apparatus. 
Obviously these results cannot replace a full Monte-Carlo simulation, for detailed studies. However they are 
exceedingly useful to improve the basic understanding, for a fast outlook and for a rough cross-check of the 
detailed Monte-Carlo simulation results. The basics of the optimisation of the experiment design can be easily 
understood using the results we will present, leaving the hard work of a full Monte-Carlo simulation to a second 
phase. In fact the very many important parameters affecting the performance do not allow to perform a full 
Monte-Carlo simulation of all the possible cases: the results presented here provide therefore valuable starting 
points for detailed Monte-Carlo simulations. 

A full Monte-Carlo simulation was developed in the framework of the EUSO Collaboration: ESAF [5J [9], 
allowing to extract detailed predictions. In fact most of the work presented in this paper is the result of the 
development of the EUSO experiment [10] , a path-finder mission of the European Space Agency (ESA) , which 
successfully completed its phase A study in 2004 but was frozen afterwards due to programmatic and financial 
reasons. A large part of the scientific community is now looking forward to both the path-finder Japan-lead 
mission JEM-EUSO [TT] and to the most challenging mission (super-EUSO [H]) in the framework of the ESA 
Cosmic Vision program 2015-2025 [IB] . 

The results presented in this paper were extensively used in the concept study of the EUSO |14] and 
super-EUSO [Q] experiments. 

In fact the design of a space-based apparatus for the observation of UHECP is a very challenging task with 
very little design margins, as it will be clear from the results we present. It is therefore important both to optimise 
the design from the beginning and to set safe design margins since the beginning of the concept study itself. 
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2 The observational approach 



John Linsley [15], in 1982, first suggested that the Earth atmosphere at night, viewed from space, constitutes a 
huge calorimeter for remotely observing UHECP (SOCRAS). Yoshiyuki Takahashi, together with John Linsley 
and Livio Scarsi, resurrected the original idea in 1995, when the concept seemed to be close to be technically 
feasible (Space Air Watch). Since then a number of proposals and studies were carried on, including the OWL [TB] 
(Orbiting Wide-angle Light-collectors) project, the TUS/KLYPVE [T7] project and EUSO [TU]. In more recent 
times the JEM-EUSO [11 Collaboration aims to propose again the EUSO concept on the International Space 
Station while the super-EUSO proposal [12] . for a most challenging next-generation experiment, has been 
recommended by ESA for technological developments. 

A space-based experiment can detect the near-UV air scintillation light isotropically produced during the 
EAS development in the atmosphere by the interaction of the EAS secondary particles with the air molecules. 
The measurement of the scintillation light allows to record the EAS development. Additional information can be 
gathered by observing the Cherenkov light diffusely reflected at the Earth surface (by land, sea or clouds). The 
Earth atmosphere plays the role of a giant calorimeter, although it is passive, continuously changing and outside 
the human control. 

This approach is complementary to ground-based observations. In fact EAS develop close to the Earth 
surface. Therefore, thanks to the large average distance from the EAS of an Earth orbiting apparatus, a large 
FoV apparatus can watch a huge target of atmosphere, so that the space-based observation is best suited for 
observing very low fluxes and/or particles weakly interacting with the atmosphere. Nevertheless, due to the 
larger distance from the EAS with respect to ground-based experiments, the signal is much fainter and therefore 
the EAS reconstruction is worse than at ground-based experiments. 

The required apparatus is an Earth- watching large aperture, large FoV, fast and highly pixelised digital 
camera for detecting near-UV single Photo-ns superimposed on a huge background, capable of three to five years 
of operation in space, at least. 

Any EAS is seen by the apparatus as a point moving on its Focal Surface (FS) with a direction and an angular 
velocity depending on the EAS primary direction with respect to the line-of-sight. These characteristics allow 
one to distinguish the EAS from the various types of background (see section [5. 10 1, because, typically, those have 
a different space-time development. 

The general scheme of the observational approach is shown in figure [l] 



3 Scientific requirements 

The following typical scientific requirements for UHECP observation from space are assumed for indication 
purposes only (see [12 ). Depending on the scientific objectives different requirements might be used by adapting 
the relations to the case considered. 

• All sky coverage. 

• Spatial granularity and resolution of the EAS image projected at the Earth: A£ ~ 0.5 km or better, in 
order to ensure a good enough reconstruction of the UHECP. 

• Sampling time of the EAS signal of order of At ~ {is or better, in order to ensure a good enough 
reconstruction of the UHECP. 

• Angular resolution on the reconstructed primary UHECP direction: Aijj < (1° ~ 3°) for a large enough 
subsample of events, in order to allow source identification and taking into account the deviation induced 
by magnetic fields on charged particles. 

• Energy resolution: AE /E < 0.3. 

• Resolution on the depth of the EAS maximum measurement A Am < 50 g/cm 2 (accounting for the intrinsic 
variability of the EAS development). 

• Energy threshold: £th ~ 1T0 19 eV, with a flat efficiency plateau at E > -Eth to keep systematic effects 
well under control. This alos ensures a fair overlap with the energy spectrum observed by ground-based 
experiments. 

• Capability to measure EAS with energies up to Emax ~ 10 21 eV. 

• Instantaneous geometrical aperture, A g , one order of magnitude larger than currently existing and/or 
planned ground-based experiments: A g > 10 6 km 2 • sr. 
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Figure 1: Air Watch observational approach: H is the orbital height, 
angle between the optical axis and the local nadir. 



7m, the FoV half-angle and otitis the tilt 
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3.1 The experimental apparatus 

The design of an apparatus for the UHECP detection from space is a very challenging task, because of the low 
expected rate of events, the faint signal from any EAS and the tight technical constraints (especially on mass, 
power, volume and telemetry) imposed on a space experiment as well as the harsh space environment. 

The engineering is very complex and the design has a very strong impact on the scientific performances. A 
careful optimisation is therefore mandatory. 



3.2 Architecture of the instrument 

The required experimental apparatus is made of the following parts. 

• The main digital camera, which is assumed to have cylindrical symmetry around the optical axis, operating 
in the near-UV, a large aperture, large FoV, fast and pixelized single Photo-n detector, consisting of: 

— the main reflective deployable optics; it consists of: 

* the main mirror: a large, lightweight, segmented, nearly spherical, deployable mirror; 

* the corrector plate on the entrance pupil which needs to be deployable as well; 

* the optical filters; 

* active control mechanism for both the mirror and the corrector plate; 

* supporting structure. 

— The Photo-Detector (PD) on the focal surface (FS) of the optics; it consists of: 

* a bidimensional array of sensors; 

* the light-collection system on the sensor; 

* the f/e electronics chip and ancillary electronics; 

* the housing; 

* the back-end, trigger and on-board data-handling electronics. 

• The Atmospheric Monitoring System (AM), capabile to monitor the relevant atmospheric properties by a 
suitable instrumentation. 

• The Monitoring, Alignment and Self-Calibration system (MAC). 

• A Radio Pulse Detection system (RPD). 

• The Central Control Unit (CCU), providing the intelligence to all the systems. 

• The system and other ancillary parts including: 

— the mechanical structure, including the external protection; 

— the thermal control system; 

— the power system, including solar panels and batteries; 

— the shutter; 

— the optical baffle. 

• The ground-support equipment. 

3.2.1 Mission requirements 

The following mission requirements are envisaged. 



An accurately designed orbit (see 5.1). The ground-track of the satellite orbit needs to be optimized 
according to the following parameters: minimization of natural and man-made light entering the FoV (Sun, 
Moon and any other light source); maximization of the rate of passages above a few fixed points at the 
Earth surface, for both exploiting ground calibration sources and hybrid observations and cross-calibrations 
with ground-based experiments. A free-flyer can give many degrees of freedom in the choice of an optimized 
orbit and the orbit strongly affects both the performance and the operations. A careful orbit optimization 
is a complex task requiring a full phase-A study. However it is not considered to be a critical issue at this 
stage as the orbital parameters can be varied to a large extent. 
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• Observation of a mass of atmosphere as large as possible, attaining a good enough geometrical acceptance 
which implies a large enough FoV. 

• No FoV obstruction and no parasitic lights coming from any other spatial device nor from the Earth. 

• Dimensions after deploying shall be dictated by the optics, including the optical baffle and the optical 
shutter. 

• Light-tightness: the interior of the instrument shall be light tight such that the parasitic lights impinging 
onto the Photo-Detector will be two orders of magnitude less than the expected night-glow background 
rate in order not to spoil the energy resoluzion. 

• Thermal control shall be provided to stabilize the large surfaces of the instrument with large power 
consumption. 

• Electrical requirements: the scientific operations will have a low duty cycle. It is envisaged that during 2/3 
of the orbit the Instrument will be in a standby status of low-power consumption. The power consumption 
quoted for operations is thus required during roughly 1/3 of the orbit only. 

• Attitude: approximate nadir pointing; pointing accuracy is not a critical factor provided the absolute 
direction of the instrument axis is known/measured for off-line use. It is assumed nadir pointing to within 
a few degrees with pointing direction known offline to within A% ~ 0.5°, well below the expected angular 
resolution of the instrument. 

• Telemetry and telecommanding: the expected event rate will be high and it will depend on the final orbit, 
affecting the event rates. A precise estimation is currently lacking. 

• Required lifetime: five years minimum, ten years goal. 

3.2.2 Requirements for the digital camera 

The basic requirements for the apparatus follow. 

• Single Photo-n detection in the wavelength range WR = {330 nm 400 ran} to detect the air scintillation 
signal. Shorter wavelengths suffer absorption from the ozone layer. 

• The faintness of the signal requires high photon collection capabilities and photon detection efficiency, 
as well as low noise (both intrinsic and externally generated). Small cross-talk and after-pulse rates are 
required to keep a good enough energy resolution. 

• A dynamic range spanning some three orders of magnitudes in the EAS energy, that is in the signal, in 
order to cover the range from the required threshold energy up to Emax- 

• An efficient and selective trigger system, to achieve a good background rejection on-board, and a powerful 
on-board data handling system. 

• Complete modularity of the apparatus to reduce the risk of single point failures. 

• Compatibility with the constraints relevant to a space mission including mass, power, volume, telemetry, 
as well as the many environmental factors. 

4 Definitions and assumptions 

We will use some simplifying assumptions in order to deduce the basic relations and use analytical or 
semi-analytical relations whenever possible. 
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4.1 Some geometry definitions 

The EAS properties depend, to a first approximation, on the primary particle identity, on its energy E and zenith 
angle, 6. The features of the EAS image also depend on the location of the EAS image inside the FoV of the 
apparatus, as defined by the field-angle 7 of a specific point of the EAS image, and the azimuth angle of the EAS 
image projected at the Earth surface. 

The angle between the latter and the radial direction in the FoV is called ip a . Let the EAS direction azimuth 
angle, with respect to a global reference system Oxyz, be ip and the azimuth angle of the line of sight in the FoV 
be 0. Clearly ip a — tp — <fi (figure [2]). The the unit vectors of the EAS direction n and of the line of sight s in the 
frame Oxyz are respectively 

h = {—sin 9 cosy, — sin # sin y, — cos#} , 
s = {sin 7 cos 4>, sin 7 sin 0, —cos 7} . 

For an apparatus with cylindrical symmetry, there is no dependence on 4>, but only on ip a . 
Let j3 be the angle between the EAS direction and line of sight s from the apparatus to the actual EAS 
location: 

cos(/3) = cos(7) cos(6>) — cos('0a) sin(7) sin(0) . (1) 




Figure 2: Definition of the angles in the FoV with respect to a global reference frame Oxyz. The point D is the 
detector position while the point P is the impact point on the EAS on the Earth surface. 

4.2 Reference conditions and general assumptions 

In this paper, unless otherwise specified, we shall use, as a reference, the conditions and parameters summarised 
111 table [U 
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Quantity 


Reference value 


Other values 


Energy 


E = 10 19 eV 




Zenith angle 


e = 50° 


9 = 30°, 70° 


Azimuth angle 


i>a = 90° 


V> a = 0°, 180° 



Orbital height of the detector H = 700 km H = (400 4- 1000) km 

FoV aperture (half-angle) 7 M = 20° 7 M = 15°, 25° 

Field Angle of observation^] 7 = 15° 7 = 10°, 20° 

Tilt angle of the detector a t ut — 0° 

Total photo-detection efficiency epD =0.1 



Operating Wavelength Range (WR) 330 nm < A < 400 nm 

Average atmospheric transmission (in WR) K atm > 0.4 

Random Background in WR [ph • m-^-^r -1 ] £ = 5-10 n B = (0.3 -=- 1.0) -10 12 



Table 1: The parameters and conditions used as reference in this chapter. 



"The field angle 7 is chosen approximatively at half of the FoV. 



The basic parameters affecting the EAS reconstruciton are: the energy, nature and direction of the EAS, the 
angle between the EAS direction and the line-of-sight from the apparatus to the EAS, (3, the field angle of the 
EAS image, 7, and the altitude of the apparatus above the Earth, H. 

Moreover the following general assumptions will be used in all the paper. 

1. The EAS energy is fixed by the scientific requirements to have a good superposition with the spectrum 
observed by ground-based experiments. Therefore a reference energy E rc f sw 1 • 10 19 eV is used. 

The EAS direction and tp a strongly affect the kinematics of the EAS image signal and therefore its 
observability and observed features. We will consider EAS with 9 « 45°, which grants a good extension of 
the EAS with a negligible shortening of the longitudinal profile due to the impact with ground. In order to 



describe an average behaviour an angle ip a — ±90° will be used as a reference (see section 5.6 1 



2. The discussion will be based on typical hadron-induced EAS, with energy E w (0.1 ~ 10) TO 20 eV and a 
simple apparatus, looking downward the Earth, whose configuration and characteristics will be defined 
during the discussion. 

3. The EAS is geometrically modeled as a point moving on a straight line at the speed of light, as, basically, 
any EAS is seen from any realistic space-detector as a unidimensional object, on all practical purposes. A 
reference EAS energy E rc { w 1-10 19 eV is used. 

For several reasons, to be discussed later, we only consider EAS with a zenith angle 30° < 9 < 70°. A very 
inclined EAS requires in any case a diferent treatment thna an almost vertcal EAS. 

One should note that when observing from space the case might arise of up-going EAS, that is EAS with 
zenith angle larger than 90°. 

4. Exponential density profile of the Earth atmosphere (assuming an isothermal atmosphere) as a function of 
the height, h, above the sea level [18] : 



p{h) = po exp 



h 
h 



with po = 1.2249 kg/m 3 and h = 8.4 km . (2) 



In some calculations we have also used, as a cross-check, the Linsley's parametrisation of the US Standard 
Atmosphere (151 EH] . 

5. The flat Earth approximation is used whenever applicable, that is when horizontal distances are negligible 
with respect to the Earth radius. When this approximation is not good enough a spherical Earth with 
radius R$ = 6371 km will be used. 

The flat Earth approximation does not substantially affect the results on the EAS development for zenith 
angle 9 < 70°. In fact, for 9 = 70°, the difference between the exact and the approximated value of the 
EAS slant depth at ground is about ~ 30 g/cm 2 , as shown in figure [3] which is smaller than the desired 
resolution on the depth of the EAS maximum: AXy 
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Figure 3: Absolute value of the difference |AX| between the true slant depth at ground and the one in flat 
Earth approximation versus the zenith angle. 



Therefore the following approximate relation between the quota above the Earth surface, h, and the distance 
t measured along a straight line with zenith angle 8 will be assumed: 



df 



1 



dh cos 9 

The exact relation for a spherical Earth, in terms of the Earth radius ii®, is: 
d£ h + R® L h 



(3) 



dh 



h 2 + 2R (f) h + i?i cos 2 i 



1 - — — tan^ y | sec t 

Ra 



(4) 



clearly showing the condition allowing to neglect the sphericity of the Earth: £, h <C Rq 

It is also worth remembering, for future uses, the d 
the appropriate boundary condition on X at some 



r / ^ dX 

It is also worth remembering, for future uses, the definition of X (the so called grammage): = p, plus 



6. The Gaisser-Hillas (GH) parametrisation [55] of hadron-induced EAS parameterises the EAS 
longitudinal profile as a gamma function: 



N(X) = N Q 



X-Xn 



X 



M 



X 



exp {(X M - X) /A] 



= iVn 



X-Xr 



x R 

for X > X Q 



Xr/\ 



exp [(X R + X Q - X) /A] , 



(5) 



where N(X), the number of charged particles at the slant depth X, is expressed in terms of No, the number 
of charged particles at the EAS maximum, Am, the depth of the EAS maximum, Xq, the starting point of 
the EAS development, Ar = Am — A and the parameter A ~ 65 g/cm 2 . The second expression is used 
to make explicit the invariance of the EAS longitudinal profile with respect to shifts of the first interaction 
point Aq. 
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The pseudo-age of an EAS at point A, is defined, closely following the age definition for an electromagnetic 
shower, as 

S[A) - (X-X ) + 2(X M -X ) ' W 

in terms of the actual slant depth and the slant depth, A , of the first interaction point. 

In our calculations (see [22]) Xq can be fixed at the first interaction point and the A parameter is rather 
energy independent and very similar for both proton and iron induced EAS in the energy range we are 
studying. Moreover Nq — aE with a ~ 0.6/GeV [53]. We can also roughly assume that for a primary 
proton EAS Xq ps 35 g/cm 2 and for a primary iron EAS Xq « 10 g/cm 2 with no significant dependence on 
the energy. 

In practice the parameter X poorly correlates with the actual first interaction point, X\. Fixing Xq = X\ 
when fitting the EAS profile with the GH parametrization, produces a worse fit with Am reduced by about 
10 g/cm 2 . 

Typical values for the EAS parameters, as taken from [55], are summarised in table [2] The values 
of Am are the average value for different interaction models (MOCCA-Internal, MOCCA-SIBYLL, 
CORSIKA-SYBILL and CORSIKA-QGSJET). Note that these values are affected by X when using the 
parametrisation of equation See [35] for details. 



Particle 


Enerj 


ry [eV] 


A M [g/cm 2 ] 


Particle 


Energy [eV] 


A M [g/cm 2 ] 


Proton 


1.0 


10 iy 


824 ± 55 


Iron 


1.0-10 19 


724 ± 20 


Proton 


5.0 


10 19 


853 ± 58 


Iron 


5.0-10 19 


754 ± 20 


Proton 


1.0 


10 20 


880 ± 54 


Iron 


1.0-10 20 


784 ± 19 


Proton 


5.0 


10 20 


907 ± 49 


Iron 


5.0- 10 20 


814 ±19 



Table 2: Typical values of EAS parameters from 22] . 

The height hyi of the EAS maximum as a function of the EAS zenith angle 9 is shown in figure [4] for three 
different values of Am using Linsley's parametrisation of the atmosphere density profile. The exponential 
profile for Am = 900 g/cm 2 is also shown for the sake of comparison: it is clear that the choice of the 
atmosphere density profile does not strongly affect hyi up to 8 « (60° 70°). 
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Figure 4: The height of the EAS maximum versus the EAS zenith angle. 
The lateral EAS profile is not accounted for, as it is hardly accessible when observing from space. 
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7. A typical hadron-induced EAS of E = (0.1 + 10) -10 20 eV with 9 < 70° is seen on the FS of an orbiting 
apparatus as a track not longer than a few degrees, as discussed in detail in section |5.6| Actually this 
depends on the instantaneous height of the satellite from ground: a few degrees length must be considered 
an upper limit, as it is calculated for H ss 300 km, which is a very small orbital height for any satellite 
requiring a reasonably long lifetime Q 

We can therefore neglect any possible change of the apparatus properties as a function of the field-angle for 
any given EAS, even if we are dealing with a large FoV optics. Obviously this would not be true for a very 
inclined EAS, which can appear many degrees long and cannot be easily described in a simple parameterised 
form. 

8. The EAS develops in the atmosphere at an height ft, < 20 km. In this part of atmosphere the 
scintillation yield in the wavelength range WR can be considered as nearly constant [24]. Roughly 
Y ~ (4.2 ± 0.2) photons particle -1 m _1 . When a better precision is required one can use the yield 
measurements of [25] . 

9. The random background (RB) is assumed to be uniform, isotropic and constant on the space-time 
scale of the EAS development. We assume for the RB in the wavelength range WR a reference value: 
B ps 0.5±^-10 12 photons • m~ 2 • s" 1 • sr" 1 . 

From the current measurements we already know that, under some circumstances, it can be up to a factor 
two larger depending on the conditions (including moon phases and cloudiness). 

It is assumed that, by continuously measuring the background nearby the EAS, both in space and in time, 
the underlying RB can be subtracted away, in real time. Due to the large RB this is an essential assumption 
in order to extract the faint EAS signal from the background at energies near -Eth- 

The possibility to subtract the background in real-time is justified by the rather large expected RB rate, 



such that the relative error on the background estimate is small (see section 6.4 1. Therefore an appropriate 
statistical estimator used to compare the signal to background is S/ \B. 

10. Atmospheric transmission is naively modeled taking only into account Rayleigh scattering. In fact Mie 
scattering is only important at low altitudes of a few km, while the EAS predominantly develops above 
a few km height (see figure [4}. Mie scattering is therefore ignored at the current level of approximation. 
Another factor which might turn out to be important is multiple scattering of the signal, which is however 
difficult to study without a full Monte-Carlo simulation. 

A most important issue is the atmospheric transmission as a function of the zenith angle direction. Table [3] 
shows the ratio, r, between the atmospheric transmission from ground to infinity at a zenith angle 9 and 
the same quantity in the vertical direction, all quantities being evaluated in the flat Earth approximation. 
The table shows that at 9 « 75° the atmospheric transmission is reduced by one order of magnitude with 
respect to the vertical transmission. This fact is a crucial one in the evaluation of the effectiveness of a 
tilting of the apparatus with respect to the nadir direction, as it will be discussed later. 

11. The atmospheric geodesic refraction is normally negligible. Indeed, according to |18j . the correction for the 
refraction when observing an object at the Earth surface from an height H, at a zenith angle x with respect 
to the nadir, is Ax « 62.37" tan x, for x ^ 80° at a wavelength A = 337 nm. 

So we have A% ~ 0.008° at \ — 25°, that is about the maximum field-angle in the case of no tilting of the 
detector optical axis with respect to the nadir. The effect of the refraction is small with respect to the pixel 
size (typically of ~ 0.1°) and it can be neglected. 

However the effect of atmospheric refraction can be huge with a tilted apparatus, when observing at large 
angles with respect to nadir. 

12. Clear sky conditions are assumed. In fact a study taking into account all the many effects of clouds requires 
without any doubt a detailed Monte-Carlo simulation. Therefore we will find best-possible results. 

13. Basic parameters of the apparatus are its altitude above the Earth, H, (or, in more general terms, the 
orbital parameters) and the tilt angle atat between the optical axis and the local vertical. 



1 This angle, as an upper limit, is easily estimated computing the EAS length between the first interaction point and the ground. 
The real angle seen by the detector is smaller due to the background predominance with respect to the signal at the beginning of the 
EAS. 
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Tabic 3: Ratio 1Z between the atmospheric transmission from ground to infinity at a zenith angle x an d the 
same quantity at \ — 0°, all quantities being evaluated in the flat Earth approximation (for a wavelength 
A = 337 nm) using equation ^ . 



The choice of the orbit height of the apparatus requires a trade-off taking into account the low signal 
produced by the lower energy EAS and the requirement to observe a target of atmosphere as large as 
possible. Moreover practical and technical constraints may limit the orbit height. 



The orbit height of the apparatus, H 7 is assumed to be: 300 km < H < 1000 km (see section 5.1 1. In fact 
orbits lower than about 300 km suffer too much drag and require therefore frequent reboosts. Orbits higher 
than about 1000 km are disfavoured by the fact that the radiation environment undergoes a substantial 
change at about 1000 km height with the inset of Van Allen belts, whose high level of trapped radiation 
may greatly influence the spacecraft. Moreover orbits higher than about 1000 km turn out to provide too 
high a threshold for any reasonably-sized apparatus. 

The orbit parameters must be tuned to optimize the expected results: energy range and instantaneous 
geometrical aperture can be tuned but constraints have to be taken into account. 

Some of the above assumptions are optimistic with respect to the real conditions of the experiment (the 
clear sky assumption, for instance). This means that the results we will obtain are necessary, but possibly not 
sufficient, requirements for observations. 

4.3 The air scintillation signal 

It is worthwhile to recall the expressions giving the air scintillation signal from the EAS. Any wavelength 
dependence is implicit in all the following formulas. 

Let to be a linear position coordinate along the EAS development (Ato = cAt), let x[w] = a;[w[i]] be the 
current EAS position and y be the current position of the observational apparatus, let Ye^M) = ^(^[w^]]) be 
the air scintillation yield in the desired wavelength range at the EAS position, let N c (x[w]) = A r c (a;[it;[i]]) be the 
number of charged particles in the EAS at w and let A be the entrance pupil area of the optics (which is actually 
the cross-section for a photon hitting the FS when 7 = 0). 

Let the atmospheric transmission from the generation point to the observation point be Tp^x^y]. The 
atmospheric absorption due to Rayleigh scattering only can be roughly parameterized as: 



T A [x[w],y] = r a [x, z] ~ exp 



Poho 



M(x) exp 



A(AW w 1 \ h 



(7) 



with ^-^ ~ 0.7 for A = 337 nm , (8) 
A 



where h w 8.4 km is the atmosphere scale height at the sea level (see the ([2|), po = 1-292 kg-m~ 3 is the 
atmospheric density at the sea level, z is the altitude of the photons emission point, A(A) is the Rayleigh mean 
free path at the wavelength A (see [55]) and M(x) is the airmass function at the zenith angle x- For x ^ 80°, 
M(x) ~ sec X; a better approximation can be found for example in |27j . 

Let the overall detection efficiency of the PD, including the f/e electronics efficiency, be £pd[7, X\\ it is is the 
probability that a photon reaching the PD will fire the f/e electronics producing a recorded photon hit. It might 
depend on many parameters, X ', and it depends on the field-angle 7, in particular. In fact at different field-angles 
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the incidence angles of the incoming photons may change leading, for instance, to a change in the intrinsic sensor 
detection efficiency. 

Let the total efficacjj^] [35] of the optical apparatus be (yleoM) and the total throughput efficiency]^] [25J be 
£o [7] ■ The previous definitions refer to all the photons reaching any point on the FS and do not include the filter 
transmission, 77f[t], which is intentionally kept as a separate parameter as the filter transmission is one the the 
key parameters to play with. 

In order to quantify the detectable signal one should build more appropriate quantities than the total efficacy 
and total throughput efficiency. In fact for any given direction, at an angle 7 with respect to the optical axis, 
only the photons reaching the FS close enough to the centroid of the distribution of all the photons focused on 
the FS are useful for the EAS reconstruction. In fact photons too far from the centroid will contribute to the 
veiling glare, that is diffused background on the FS. This effect might be particularly dangerous because, thanks 
to the large FoV, it is most likely that bright sources are present somewhere in the FoV: in case of a large veiling 
glare these bright sources might significantly increase the background level on the whole FS. 

A suitable fiducial region (called f2), around the centroid of the distribution on the FS of all the photons, 
must be defined. The fiducial region £1 may be defined in different ways. The photons on the FS will then be 
classified into three classes, as follows: inside the fiducial region f2 (transmitted photons); outside the fiducial 
region Q but detected somewhere on the FS (photons giving rise to the veiling glare); any other photon (lost or 
absorbed somewhere). It is appropriate to adopt as a fiducial region a region of the order of the spot-size and/or 
pixel size. 

Let eb[t] be the fraction, out of all the photons on the FS, falling inside the desired bucket size: this is 
actually a measure of the encircled energy fraction. 

Define then the triggering efficacy as Aeo[7]£b[7] = -Ae^M an d the triggering throughput efficiency as 
£o[j}£Bh] = e'ob]- 

Let AS[w] (with the subscript f2) remind that the detected signal depends on the region on the FS around 

the the centroid of the spot (bucket) where photons are considered to be detected photons instead of background. 
The photons outside just contribute to the veiling glare. 

The number of signal photons produced by a segment Aw of the EAS, at coordinate w, and detected by the 
apparatus is thus: 

AS[w] = e \ l \ C 1 J 2 T A [x, y] (AeoWs^M^epuh X] . (9) 
n V 47r (x[w\ — y) J 

AS\w\ with the subscript f2 is meant to remind that those equations are only meaningful after specification 

of the region on the FS over which photons are considered to be detected photons, that is photons that are 
close enough to the centroid of the image spot. 

It is worth noting that when looking at faint signals (a few photons per pixel) the errors in formula [9] will be 
typically dominated by the Poisson statistical fluctuations, at least for the faintest EAS. In fact it will be shown 
in section |6.3| that for the faintest EAS, with a number of detected photons N w 100, one has less than ten 
photons per pixel with a relative statistical error of the order of 0.3, larger than the expected systematic errors 
(~15%)0 

Using the simplifying assumptions listed in section [4] it is easy to build a parameterised semi-analytical 
description of the relation [9] which will be used in the rest of the paper (cfr [251 13"U]V 



AS[i 



A EP N c e T 



D 
Do 



YAw , (10) 



where ex is the total efficiency (optics plus photo-detector and filter), D the distance from the detector to the 
fluorescence emission point (it can be assumed as a constant) and Dq ~ 12 km is the extinction length of light 
due to the atmospheric scattering (at A = 337 nm). 



4.3.1 Single-photon detection factors 

The atmospheric transmission Ta[^M, y] is the first basic single photon detection factor. It is out of control but 
the experiment design should be optimised to minimise its effects. 

2 It is defined as the ratio the number of photons incident on the entrance pupil per time unit with an angle 7 and the photon 
irradiance [ph-m — 2 s -1 ]. 

3 It is defined as the ratio the number of photons incident on the entrance pupil per time unit with an angle 7 and the total photon 
flux collected by the pupil (that is given by the photon irradiance times the entrance pupil area). 

4 Obviously when fitting the profile with a given functional dependence the resulting error will be reduced by the fitting procedure. 
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The photon collecting power of the optics is: A£o[t]£b[7]?7f[7]> including the effect of PSF and filters, ft 
gives the conversion from the incoming number of photons per squared meter, at an angle 7 with the optical 
axis, to the number of photons reaching the FS inside the region £1. It is the quality factor for the optics as far 
as the number of collected photons is concerned. It is the basic and unique sizable parameter to define the optics 
photon collecting power. 

The overall detection efficiency of the PD, £pd[7j X], is another basic parameter which depends on a large 
number of factors. These factors are often very close to one but, when multiplying a large number of them, the final 
result may significantly depart from one. The Photon-Detection-Efficiency of the sensor is the most significant 
parameter which is typically much smaller than one. For standard PMT, which have already been tested in the 
space environment and proved suitable for applications in space, £pd — 0.12 4- 0.15. Recently several new high 
quantum efficiency photo-detectors became available on the market: GaAsP Hybrid PhotoDiodes (HPD) |3T] and 
Geiger- Avalanche Photo-Diode (GAPD) j32j[33]. These devices, exploiting the intrinsic quantum efficiency of 
the solid state device, aim to reach a £ppj = 0.4 -j- 0.6. 

The overall capability to detect photons of the instrumental apparatus can be expressed through the 
photo-detection efficacy, defined as 

£p D = A ■ e' [7] • tjfM • epdH ■ (H) 
4.3.2 Multi-photon detection factors 

Other more complex factors will affect the EAS detection, involving multi-photon correlations and correlations 
between neighbouring pixels, such as: the read-out electronics and trigger efficiency, £tpj, and the analysis and 
event reconstruction efficiency, £ana- These effects are very difficult to estimate without a full Monte-Carlo 
simulation j8j [9j [34j [35] , as they depend on the contribution by all the photons at once. 

Therefore naive expectations coming from single photons studies can be far too optimistic: once more simple 
estimates may turn out to be optimistic. 

5 Design and optimisation of a space-borne apparatus 
5.1 Orbit 

The orbit parameters can be tuned to optimise the expected performances, in particular the energy range and 
the instantaneous geometrical aperture. 

The orbital height (semi-major axis of the ellipse) is one of the most important parameters. In fact an higher 
altitude implies a larger mass of observed target atmosphere (and therefore a better instantaneous geometrical 
aperture) but also an higher energy threshold because of the smaller EAS signal. 

Varying the orbital height is useful in order to span a larger range of energies: with an elliptic orbit the 
satellite spends more time at a higher altitude, gaining in effective aperture, but it also spends some time at a 
lower altitude, decreasing the energy threshold, where a long measuring time might be not necessary. 

An orbit height spanning the range 400 km w Hmin S H < Hmax ~ 1000 km, for instance, would give 
the capability to extend the energy range, with respect to a fixed height, by taking data at different heights, 
that is shifting up and down the UHECP energy range with a scaling factor: (^max/^min) 2 ~ 6. Moreover 
the instantaneous geometrical aperture at apogee would be the same factor larger than at the perigee. On the 
positive side one should note that the angular granularity, at fixed AL, scales as H~ l , and not as H~ 2 . 

Another way to vary the height might be to use different almost circular orbits during the mission lifetime 
(for example part of the time at a lower altitude and part of the time at an higher altitude) . Natural orbit decay 
might be exploited as well. 

One should remember that the orbit lifetime depends on the height and it strongly depends on the epoch of 
the Solar cycle (that is on the epoch of the Mission) , which affects the atmospheric density. 

Another key point to keep in mind is that the ballistic coefficient of the satellite is expected to be rather low, 
due to the large expected area-to-mass ratio. Therefore, very roughly, the orbit lifetim^Jis expected to be of the 
order of a few days for H ps 300 km, a few weeks for H ps 400 km and some years at H sa 700 km (see figure [5]) . 

The orbit inclination should be chosen in order to span as much as possible the Earth surface in order to 
obtain a sky coverage as uniform as possible. 

An important parameter that depends on the orbit is the duty cycle, that should be as large as possible. For 
this purpose, the orbit design must be further optimised in order to avoid as much as possible the light pollution 

5 Formula and data for calculation can be found in 36 . 
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Figure 5: Lifetime L (in months) of a satellite (mass ~ 3000 kg; total area ~ 150 m 2 ) versus the orbital height 
H. Solid line: solar cycle minimum; dotted line: solar cycle maximum. 

from the man-made sources at ground and from the lightings (that are concentrated mainly on lands) and the 
auroras. 

Note also that the effects of a short duty cycle can be recovered increasing the mission lifetime. 

The orbit design can be also optimised in order to access as much time as possible over fixed targets at ground 
(ground EAS experiments, weather stations, ground-based calibration sources, etc.). In particular, the accesses 
to other EAS experiments, namely the Pierre Auger Observatory, are very useful for the apparatus calibration. 

5.2 Optical system design 

The two basic optical parameters, affecting both the performance and the engineering of the experiment, are the 
Entrance Pupil (EP) area, Aep, and the Field of View (FoV) of the optics, defined as the half-angle of the FoV 
circular cone, 7m- 

The large desired aperture calls for an //# as small as possible in order to reduce the size of the FS. However 
the large desired aperture implies, in any case, a large FS size. 

5.2.1 Entrance pupil of the optics 

The apparatus sensitivity to low signals is affected by a basic and unique sizeable parameter: the optics aperture, 
that is the entrance pupil area, A EP . Increasing the aperture, to within external constraints, gives a guaranteed 
improvement of the performance: it is just a matter of technology and money to build such a large apparatus. 

As a very large aperture optics is needed a deployable optics is most probably required. 

The optics aperture, also fixes the size of the apparatus. The lower limit for A^p is basically set by the 
requirement on the energy threshold, while the upper limit is fixed by external constraints. 

The total number of detected photons is roughly proportional to solid angle subtended by the optics entrance 
pupil as it is seen by the EAS: w A E p/H 2 . Therefore, for a required energy threshold, the choice of A^p and 
orbit height is driven by the number of photons that one wants to observe. 
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5.2.2 FoV of the optics 

The large desired aperture implies in any case a large FS size. 

In order to gain instantaneous geometrical aperture, either a larger orbit height or a tilt with respect to nadir 
might be more effective than increasing the FoV. In fact these options might have two beneficial effects: to keep 
the FS size limited and to ensure a better optics performance, as the performace of the optics tends to become 
worse, normally, when the FoV increases. In fact the large required FoV and the necessity to keep the optics 
simple typically give a PSF which is not good enough to detect the lateral EAS dimensions. Even if one manages 
to get a good enough PSF the required number of channels would increase unacceptably. In our discussion we 
will consider a FoV aperture (half-angle) 15° < 7m < 25°. 

5.2.3 Reflective or refractive ? 

There are virtues and defects for using refractive versus reflective optics in this kind of experiment. A refractive 
system is possibly difficult to deploy and a complex supporting structure is required due to the large dimensions. 
A reflective system would suffer from the obscuration due to the FS and the very limited room for the FS 
Photo-Detector. The latter might be a real issue, depending on the kind of sensors one is going to use. 

Ona main drawback of a reflective system is that the curvature of the FS (that is center of curvature of the 
FS on the opposite side with respect to the incoming light) is such that the filling of the FS with the array of 
sensors is more difficult and the filling factor tends to be worse than for a refractive system, where the center of 
curvature of the FS is on the same side with respect to the incoming light. 

5.3 Effect of the observation angle (location on the FS) 

At the edge of FoV (7 ~ 7m) the EAS triggering and reconstruction is much more difficult than near the center 
of the FoV. In fact due to the larger EAS distance the irradiance (energy/photons per unit area per unit time) 
decreases as » cos 2 7. Moreover the effective entrance pupil typically decreases as cos 7: this is the so-called 
obliquity factor of the optics, due to the fact that the size of the entrance pupil is seen reduced by a factor cos 7. 
As a result the number of photons received on FS by an EAS of fixed energy scales at least as cos 3 7. 

In general the real optics triggering efficacy/efficiency, EoMsbM, also decreases while increasing the 
field-angle. This also includes the effect of the absorption losses inside the refractive materials. 

Moreover an EAS detected at the edge of the FoV would suffer from a larger atmospheric attenuation Ta [7] 
(see equation ([7])) due to the longer and more inclined path. 

As a smaller number of photons is detected from any EAS when it is observed at the edge of FoV, with respect 
to the same EAS observed at the center of FoV, one can roughly say that the FoV center (more photons from 
the EAS, smaller geometrical acceptance) is better for low energy events and low photon flux while the FoV edge 
(less photons from the EAS, larger geometrical acceptance) is better for high energy events and high photon flux 
rates. 

5.3.1 Signal and RB roll-off with field-angle 

It is important to compare the dependence of both signal and RB on the field-angle. 

Let us assume an ideal optics, that is an optics such that all the photons hitting the entrance of the optical 
system are focused inside the desired bucket. For such an optics: EqH = £ o[7] £ b[7] = cos 7. Beware that for a 
real optics £^[7] might typically drop faster with 7. 

The field-angle dependence for a given source is given by 

S[i\ ~T A [7] e ' [7] cos 2 7 ^ 7 (0 ^25°) : (1.0 -> 0.72) . (12) 

On the other hand the field-angle dependence for the RB is given by 

B fr] ~ e 'oh\ — ► 7(0 -> 25°) : (1.0 -► 0.87) . (13) 

Both signal and RB decrease at increasing field-angles, but in a different way. In fact the field-angle dependence 
for the signal to RB ratio is: 

%} ~T A [7] cos 2 7 — > 7 (0 ^25°) : (1.0 -> 0.79) . (14) 

The dependence on the field-angle of both the signal and the signal-to-background ratio renders useless the 
increase of the optics FoV above a certain extent for a nadir pointing instrument. Similar considerations apply 
to tilted instruments. 
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Many parameters (such as the trigger settings) require tuning as a function of the field-angle, that is as a 
function of the radial distance on the FS. 

One important result is this that, as the efficiency curve as a function of the energy roughly scales in energy as 
the inverse number of signal photons detected, a detection energy threshold at least sa 2 times higher is expected 
for events detected at 7 = 25° with respect to events detected on-axis [5H [35] . 

As a conclusion the performances drop significantly at large angles from nadir due to the larger EAS distance, 
the larger atmospheric absorption and the worse optics photon collection efficiency (the latter actually depends 
on the angle from the optical axis, not from nadir). 



5.3.2 Other considerations 



The FoV, together with the satellite altitude, determines the geometrical aperture (discussed in the section 5.7 1 
and therefore the number of events to be detected. The Z?ep defines the detector collection area and its sensitivity, 
but also fixes the apparatus size and, consequently, the mass and the volume of the payload. 

The diameter Z?pd of the FS is a function of the FoV (half-angle 7m) and of the focal length F: 

D PD = 2Fsm7M (15) 

If we introduce the optics /-number, //# = F/D^p then 

£pd =2//#£ E psin 7M (16) 

Note that for 7m = 25° and //# = 1 the FS is as large as the entrance pupil. A further discussion of these 
arguments can be found in |37j . 

As the apparatus is basically photon-limited, the entrance pupil shall be as large as possible, up to the limit 
allowed by external, technological and practical constraints. Therefore the size of the optics will be approximately 
set to the maximum allowable size. Therefore a reasonable assumption for any real apparatus is that the FS is 
not larger than the optics: if the latter is approximated by the optics entrance pupil, as it would be desirable for 
a high-efficiency optics, this implies the rough estimate: 2//#sin7M < 1 



5.4 The photo-detector on the FS 

As a consequence of the limits on /-number and on the pixel size, a very large area FS is required: the choice of 
the sensor has to account for this. This basically means that compact and light sensors are required. 

The pixel size is driven by two competing requirements: going far away from the Earth (or tilting) requires a 
smaller pixel size and increasing the optics aperture (and therefore the focal length) implies a larger pixel size. 

The overall photo-detection efficiency of the photo-detector is one of the most important parameters affecting 
the performance. However, needless to say, it cannot increase larger than one. Typical, realistic overall 
photo-detection efficiencies in the range (0.1 -j- 0.2) can be quoted for these kind of apparatus. Even in the 
most optimistic approach it is not realistic to expect more than a factor sa 3 improvement, in the near future, 
provided huge efforts can be devoted to the development of suitable sensors with higher Quantum Efficiency. In 
fact the main source of inefficiency, as of today, is typically the rather low Quantum Efficiency, affecting the Total 
Photo-Detection efficiency of the sensors. However also the geometrical acceptance and filling factor of the arrya 
are a crucial issue for this kind of application. 

A lot of other factors affects the overall efficiency. They are however all factors already very close to one 
(typically sa 0.9), but the product of many of them can make the overall efficiency drop. In order to make a 
significant change all of them should be substantially improved because they are already very close to one: there 
is little hope to gain anything substantial. 



5.4.1 The parameters of the Photo-Detector 

Approximate relations, useful to estimate the relevant photo-detector parameters, are summarised in this section. 
The relations in this section will be essentially geometrical, more refined estimates would require a precise design 
and full simulations. 

The characteristics of the optics have a strong impact on the photo-detector design. The Point Spread Function 
(PSF) of the optics has to match, approximately, the photo-detector pixel size. In fact a finer granularity would 
allow a better reconstruction, provided enough photons are collected, possibly allowing to measure the PSF itself. 
However a tradeoff with cost and complexity (driven by the number of channels) is unavoidable, on the other 
had, obviously, a pixel size much larger than the PSF would not exploit all the optics performance, wasting the 
efforts for the optics. 
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The FS shape and dimensions are determined by the optics. At this stage, therefore, only rough estimates of 
the corresponding photo-detector parameters will be attempted. 

A more careful evaluation can only be carried on once the design of the optics can be better defined. 

The desired number of pixels of the photo-detector can be estimated, from the desired parameters, by the 
relation 

nH 2 tan 2 7 

N S3 ^— . 17 

AL y ' 

The required photo-detector pixel size, 5, corresponding to observing a length AL on the Earth surface, can 
be estimated by the relation 

f AL 

<s«^r ' (18) 

where / is the focal length of the optics. 

The photo-detector surface has to approximate the focal surface of the optics. The latter can be assumed, 
to a first approximation, to have a spherical shape with radius equal to the focal length of the optics, /, and 
maximum angular aperture f3. Its area is then given by the relation 

v4 dct = 2 7 r/ 2 [l-cos/3] . (19) 

The FS maximum diameter is given by 

D f = 2f sin P . (20) 

Note that, at least to a first approximation, (3 ~ 7m- 

The approximate maximum number of pixels which can be fitted on the FS of the photo-detector is given by 
the relation 

Alternatively the desired pixel dimension on the FS can be estimated, in terms of the FS parameters and the 
desired number of pixels one wants to see at the Earth, by the relation 



<~i¥ ■ (22) 

Note that, to the present level of approximation and with the present parameters, the two relations |18| and 
[22] are roughly equivalent, given the relations |P7| and |19} 

The optics also determines the distribution of incidence angles of the photons on the FS, which has some 
impact on the photo-detector design. The marginal ray angle is determined by the /# of the optics and is given 
by the relation 

tan 9 max S3 — — . (23) 
2/# 

The angular granularity of the photo-detector, Aa, is given by the relations 

Aa S3 ~ - . (24) 

The solid angle coverage of every pixel, Af2, is given by 

AL 2 S 2 , 
Aft*— ^-^[Aaf . (25) 

The approximate defocusing in the direction parallel to the FS, Aw, produced by a small displacement Az in 
the direction perpendicular to the FS, is given by the relation 

Az 

Aw S3 Az tan max = — — < S . (26) 
2/# 

In order to reduce the effect of defocusing on the large FS a good fit between the ideal optics FS and the 
real FS must be implemented. In general it is not trivial to accomplish this because the optics FS has a complex 
geometrical shape. 
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5.5 Number of detected photons, efficacy and energy resolution 



The required energy resolution (AE/E ~ (20% 4- 30%)) calls for a relative error due to the Poisson statistics 
of the detected photons not larger than oat /A" ~ (0.10 ~ 0.15), assuming an equal contribution from statistical 
and systcmatics errors, as it is appropriate for a well-designed apparatus. This implies that at least one hundred 
photons must be detected from any EAS at the lowest energies, with full triggering efficiency. Obviously the 
number of detected photons impacts significantly on the other observables as well. 

Assuming an ideal optics and the reference EAS one can easily derive the requirements on the optical triggering 
efficacy from the requirement on the energy resolution. As the angular extension of the reference EAS is short 
on the FS (see section 5.6 1 the properties of the apparatus can be considered uniform. 

The time-integrated irradiance of the signal reaching the instrument, as a function of 9, is shown in the 
figures [6] and [7] From the time-integrated irradiance we can infer the required overall triggering efficacy and 
the minimum apparatus diameter (sec section [6]) . This result does not depend on the implementation of the 
optical system and the photo-sensor but some assumptions are required in order to proceed. The required overall 
triggering efficacy for observing the same EAS obviously scales as H 2 . 




5.6 The length and duration of the visible EAS image 

The angle subtended by the visible EAS image on the FS of the apparatus as well as its time duration, which are 
fundamental parameters in the experiment design, are easily estimated by determining the first and last detected 
points of the EAS. 

If one assumes to be able to subtract the background, the shape of the detected photon hits is well described 
by the Gaisser-Hillas function ^ , which is a Gamma distribution in the variable X = X — Xq : 

9(X) = f^j^X*- 1 exp (-X/p) cx A< x */ A > exp (-X/a) (27) 

where a = + 1, /3 = A, X R = X M - X . (28) 
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Figure 7: Signal time-integrated irradiance of an EAS as a function of 9 (H = 700 km and ip az = 90°). Solid line: 
7 = 15°; dotted line: 7 = 20°; dashed line: 7 = 10°. 



The mean and the standard deviation of this gamma distribution are therefore 

(X) = a/3 = X M -X Q + X 
crx = PVa = y/\(X M -X q + X) 
As we r equire to observe at least N ~ 100 photons in order to get a good enough energy resolution (see 



section 



5.5), the minimum and maximum value and the range of X in the distribution are: 

min (x) w (x) - 2er; 



max ( X ) » (x) + 3a 



range 



\--/ —x 
\ / +3(7 x 

X ' 



(Xj « 5cr 



as it can be quickly determined via simulations of gamma distributions using the physical parameters in 
section 14.21 

It should be noted that for N = 100 the range estimation is 5<7 for a gaussian distribution also but in that 
case the distribution is obviously symmetrical with respect to the mean value. However the gamma distribution 
arising from the current physical parameters has a relatively small skewnes^] and kurtosis exces^] and therefore 
it is very similar to a gaussian. 

Other values are summarized in table [3] The results show that, as expected, the range of the observed 
photon distribution (that is the observed image length) does not change by a huge amount when changing the 
photon detection capability. The dependence on the number of observed photons can be fitted by a second degree 
polynomial. 

Simulations with different parameters show that, as expected, the range of the observed photon distribution 
(that is the observed image length) is not strongly affected by little changes of the physical parameters of the 
EAS. 



6 The skewness of a gamma distribution is given by 71 = 2/y/a 0.55 for a 1-10 19 eV proton EAS. 
7 The kurtosis excess of a gamma distribution is given by 72 = 6/0 ai 0.46 for a MO 19 eV proton EAS. 
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(3.05±0.01)cr 
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N = 


= 10000 


(7.83±0.01)cr 


2.7a 


5.1a 



Table 4: Range estimation for a gamma function with parameters corresponding to a typical UHECP induced 
EAS. 



One must also take into account that a nearly vertical EAS is truncated by the ground. 

It is now easy to estimate the angle £eas subtended by the EAS image on the FS and the EAS image duration 
Teas on the FS. The dependence on the EAS zenith angle 9 are shown respectively in figures |8|f|23| 




Figure 8: Angular extension of an EAS on the FS as a function of 9 (H = 700 km and 7 = 15°). Solid line: 
^ az = 0°; dotted line: Vaz = 90°; dashed line: Vaz = 180°. 



As it is clear from these figures, £eas ^ 4°, justifying the assumptions of section [4] Note also that, for 
9 < 50°, the EAS are truncated because the EAS hits the ground, as it is evident from the cusps in the figures. 
Nevertheless the truncation does not affect the reference EAS, for which we requested N = 100 photons. Some 
numerical values for the reference EAS are given in section |6.1| 

It is obvious that the time duration of the EAS image is almost independent on both H and 7 for a fixed 

cos7/-ff 



EAS geometry (figures 11 
(figures [20| and [22|. 



13 21 and 23 1 while the angular length, for the same EAS, will scale as 



5.7 Aperture 

5.7.1 Area observed at the Earth 

In the case of a nadir pointing apparatus the geometrical area spanned by the FoV at the surface of the Earth is 
the area of a spherical cap, 

A = 27rE|(l-cos/3 M ) , (29) 
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Figure 9: Time duration of an EAS on the FS as a function of 9 (H = 700 km and 7 = 15°). Solid line: ip az = 0°; 
dotted line: ip az = 90°; dashed line: tp az = 180°. 




6 [deg] 

Figure 10: Angular extension of an EAS on the FS as a function of 9 (ip az = 90° and 7 = 15°). Solid line: 
H = 700 km; dotted line: H = 400 km; dashed line: H = 1000 km. 
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Figure 11: Time duration of an EAS on the FS as a function of 9 (V> az = 90° and 7 = 15°). Solid line: H = 700 km; 
dotted line: H = 400 km; dashed line: H = 1000 km. 
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Figure 12: Angular extension of an EAS on the FS as a function of 9 (H = 700 km and V'az = 90°). Solid line: 
7 = 15°; dotted line: 7 = 20°; dashed line: 7 = 10°. 
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Figure 13: Time duration of an EAS on the FS as a function of 8 (H = 700 km and tp az = 90°). Solid line: 
7 = 15°; dotted line: 7 = 20°; dashed line: 7 = 10°. 
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Figure 14: Angular extension of an EAS on the FS as a function of ip az (H = 700 km and 7 = 15°). Solid line: 
9 = 50°; dotted line: 9 = 30°; dashed line: 9 = 70°. 
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Figure 15: Time duration of an EAS on the FS as a function of -0az (H = 700 km and 7 = 15°). Solid line: 
6 = 50°; dotted line: 9 = 30°; dashed line: 9 = 70°. 
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Figure 16: Angular extension of an EAS on the FS as a function of ip az (8 = 50° and 7 = 15°). Solid line: 
H = 700 km; dotted line: H = 400 km; dashed line: H = 1000 km. 
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Figure 17: Time duration of an EAS on the FS as a function of ?/>az (0 = 50° and 7 = 15°). Solid line: H = 700 km; 
dotted line: H = 400 km; dashed line: if = 1000 km. 
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Figure 18: Angular extension of an EAS on the FS as a function of ip az (if = 700 km and = 50°). Solid line: 
7 = 15°; dotted line: 7 = 20°; dashed line: 7 = 10°. 
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Figure 19: Time duration of an EAS on the FS as a function of V'az (H = 700 km and 9 = 50°). Solid line: 
7 = 15°; dotted line: 7 = 20°; dashed line: 7 = 10°. 
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Figure 21: Time duration of an EAS on the FS as a function of H (tp az = 90° and 7 = 15°). Solid line: 9 = 50°; 
dotted line: 9 = 30°; dashed line: 9 = 70°. 
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Figure 22: Angular extension of an EAS on the FS as a function of 7 (H = 700 km and ip az = 90°). Solid line: 
7 = 15°; dotted line: 7 = 20°; dashed line: 7 = 10°. 
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Figure 23: Time duration of an EAS on the FS as a function of 7 (H — 700 km and i/'az = 90°). Solid line: 
7 = 15°; dotted line: 7 = 20°; dashed line: 7 = 10°. 



in terms of H, the orbital height, R®, the Earth radius, 7m the half-angle FoV of the apparatus, and /3m, the 
angle at Earth center between nadir and the FoV border which is given by: 



. . i?© + H . 
(i M = arcsm [ — sin 7 M 



7m 



(30) 



From A it is easy to estimate the total target mass, M, using the value of the vertical column density of air 
(~ 1033 g/cm 2 ). In tabled some numbers are given for different values of H and 7m- An approximated value for 
the area, using the flat Earth approximation is 



Ao = irH tan 7m 



(31) 



H [km] 


7M [° ] 


A [M0 5 km 2 ] 


M [1-10 15 kg] 


400 


20 


0.67 


0.7 


700 


20 


2.07 


2.1 


1000 


20 


4.25 


4.4 


700 


15 


1.11 


1.2 


700 


25 


3.43 


3.5 



Table 5: Area observed at the Earth and atmosphere mass target. 



5.7.2 Instantaneous geometrical aperture 

The instantaneous geometrical aperture is defined as: 



0(0, if) -ndftdA 



AJQ 



(32) 



in terms of the normalized velocity vector, z)(0, ip), of the EAS and the normal unit vector to the surface, h. 



29 



The effective aperture, for a fiat Earth, of a nadir pointing apparatus, only considering EAS which reach 
ground inside the FoV, is approximately given by the relation: 

A c = i] Q r] c (1 - T doad ) A g « i] i] c (1 - T dcad ) tt 2 H 2 tan 2 7 M , (33) 

in terms of the orbital height, H, the half-angle FoV 7m, the observational duty cycle rj , the dead time r dead 
and the cloud coverage efficiency (r] c ~ 0.5), quantifying the effect of real cloud coverage on the EAS detection 
efficiency. The expression above also gives the effective high energy (asymptotic) aperture, that is the effective 
aperture when the total detection efficiency equals one. 



5.7.3 Tilting of the apparatus 

The instantaneous geometrical aperture can increase if the apparatus is tilted with respect to the local nadir 
by some angle amt- In order to evaluate it one needs to estimate the area observed at ground with a tilted 
apparatus. As the geometrical solution of the problem to find the intersection area of the circular FoV cone with 
the spherical Earth surface is not trivial at all, the easiest way is to use a simple Monte-Carlo integration. 

The Monte-Carlo results for the intersection 3X6 el ; elS el function of a t iit , are shown in the figures [24] and [25] 

It should be noted that the corresponding horizon angle is 

{70° at H = 400 km, 
64° at H = 700 km, , 
60° at H = 1000 km 

which must be taken into account when considering the possible titling angles. 




Figure 24: Area seen by the detector at the Earth surface versus atut for 7m = 20°. Solid line: H = 400 km 
(Aior ~ 70°), dotted line H = 700 km (/3 hor w 64°), dashed line: H = 1000 km (/3 hor « 60°). 



Tilted mode allows to tune the instantaneous geometrical aperture, up to a factor (3-^5). While titling the 
apparatus is an effective way to increase the instantaneous geometrical aperture this is not the only important 
issue which comes into play. In fact the main drawbacks of tilting are that one is observing EAS at a larger and 
larger distance when looking at the far extreme of the FoV, with a drastically increasing atmospheric absorption. 
This means that the effective energy threshold in the far part of the FoV increases. Moreover the FoV is highly 
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Figure 25: Area seen by the detector at the Earth surface versus atut for H — 700 km. Solid line: 7m = 15°, 
dotted line 7m — 20°, dashed line: 7 M = 25°. 



non uniform and the angular resolution at the far extreme becomes worse unless the pixel size is reduced. The 
large FoV of the optics, together with the large atmospheric target observed, would require an excellent stray-light 
control for a tilted apparatus due to the large amount of light entering the FoV. Tilting, together with the large 
FoV, would also strongly affect the duty cycle, as the large area observed at the Earth would more often include 
day-time areas. 

On the other hand increasing the orbit height seems to be a more effective method to increase the instantaneous 
geometrical aperture while minimizing the drawbacks mentioned above. In fact one would have no losses from 
the increased atmospheric transmission and a more uniform FoV and would benefit from smaller required FoV 
in the optical design. 



5.7.4 Duty cycle 

It is not easy to estimate the duty-cycle of a space-based apparatus. A dedicated measurement is possibly 
required. 

The major limitation to the apparatus observational duty cycl^] comes from the Sun light and the Moon 
light. The fraction of time in which the detector is unable to operate due to the Sun or the Moon depends on 
the orbital elements and, during the moontime, on the maximum background rate the detector for which data 
taking is still possible. For an orbital height H ~ 400 km, an orbital inclination of ~ 50° and requiring both the 
Sun and the Moon to be safely below the horizon, the average duty cycle is r/ ~ 13%. If we accept an additional 
background flux of 100 ph m -2 ns _1 sr _1 due to the moonlight, we have rj ~ 19% [55] . 



The duty cycle may be influenced by other man-made or natural sources. They will be discussed in section 5.10 
Observation duty cycle is driven by the acceptable background level. Therefore it is, in principle, energy 
dependent but a very good knowledge of the apparatus is required. 



5.8 Pixel size and angular resolution 

Any EAS will be seen as a point moving inside the FoV with a kinematics (direction and an angular velocity) 
determined by the EAS direction relative to the line of sight from the instrument to the EAS instantaneous 

8 The observational duty cycle is defined as the fraction of time with the apparatus on, open and taking data. 
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position. The direction of the EAS velocity vector, as seen by the detector, can be decomposed into two 
components: one parallel to the line of sight and the other lying in the plane perpendicular to it. The EAS 
develops approximately at the speed of light and its distance can be considered as a known value, in the case of a 
space-epxeriment . Therefore both of them can be reconstructed from the two-dimensional image on the FS plus 
the timing information. The latter is inferred from the direction and angular velocity of the EAS inside the FoV. 

The pixel size, 6, is driven by the Scientific Requirements and constrained by available resources. In particular 
it affects: 

• the EAS reconstruction efficiency and the signal contamination: the number of signal photons divided by 
the number of RB photons on the pixel roughly scales as S/B ~ 1/8 for a pixel much larger than the EAS 
track width on the FS (which basically depends on the optics PSF) while it saturates to a constant for a 
small enough pixel size; 

• the angular resolution; 

• the Xm resolution. 

It should be kept in mind that due to the relatively small EAS transverse dimensions the EAS image transverse 
dimensions on the FS will be determined by the PSF only. 

A pixel size much larger than the optics PSF would spoil the angular resolution. On the other hand a pixel 
size smaller than the optics PSF might be useful, in case one has many photons per pixel to improve, via a 
suitable fitting procedure, the angular resolution. 

A pixel size roughly of the same size of the optics PSF turns out to be, usually, a good compromise. 

The number of pixels has a strong impact on the Instrument budgets and complexity. A trade-off on the pixel 
size is therefore very important. 

An approximate and simplified analysis, leading to determine the required pixel size, is presented below. Note 
that the following elementary analysis ignores the effect of the background (making the angular resolution worse) 
and assumes to have an unbiased estimator of the EAS arrival direction. Therefore the requirements which 
are derived must be considered as necessary requirements. On the other hand the use of the diffusely reflected 
Cherenkov flash might improve the angular resolution. 



5.8.1 Angular resolution perpendicular to the line of sight 

The expected angular resolution A/3j_ on the EAS direction perpendicularly to the line of sight is readily estimated 
by assuming to perform a linear fit. The error on the angle can be calculated from the standard relations for a 
linear least squares fit as: 

AP ± = ^= — -L= , (34) 
where N is the number of detected photons, is the standard deviation of the observed longitudinal photon 



distribution along the EAS image (easily determined from the methods presented in section 5.6 1 and 6 is the 
uncertainty in the position on the FS which can be roughly taken as the pixel size. 

Using again the sample EAS (table [2]) and exploiting the Gamma-like shape (Gaisser-Hillas function) of the 
EAS longitudinal profile, the relation between the observed (angular) EAS length, £eas> (i-e. the range of the 



sampled values) and turns out to be £eas ~ 5<j£ (for N — 100), as discussed in section 5.6 Therefore one 
finds that, in order to reach an angular resolution of the order of Aj3± ~ 1°, a FoV granularity of 6 ~ 0.1° is 
required. 

Note that the above result is consistent with the naive estimate: 

a ^~(Lts' ■ (35) 

The angular resolution A/3j_ as a function of 6, is shown in figure |26) 



5.8.2 Angular resolution parallel to the line of sight 

The relation between the observed angular velocity ojsas an d the angle f3 between the EAS velocity vector and 
the line of sight is the well known relation [501 

^eas = ^ ( 1 „™ h / ) = ^ tan [ - ) (36) 



D V sin/3 I D V 2 
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e [deg] 

Figure 26: Angular resolution A/3± as a function of 8 in the reference conditions. Solid line: H = 700 km; dotted 
line: H = 400 km (the line for H — 1000 km is very close to the first one). 



where c is the speed of light and D is the distance of the EAS. 

In the present case, one assumes that the EAS develops in the lower layers of the atmosphere, within ~ 15 km 
from the ground, so that D is approximatively known (the relative error is AD/D < 0.05). Due the non-linear 
relations between u>£j^g and (3 the estimation of the best fit error is more complex in this case. Therefore the 
error on (3 is estimated by assuming a simple measurement of the angular velocity of the EAS. One obtains: 

A/^^W^ + ^Y (37) 

WEAS \?EAS 1 EAS/ 

For the typical EAS of table [2] by assuming an EAS sampling time not larger than St ~ 2.5 fjs, the second 
term is smaller than the first one. One finds A/3y » 0.15 rad « 10° by assuming again S ~ 0.1°. One might assume 
that, with a best fit, this results will scale as TV -1 / 2 obtaining the desired A/3y 0.015 rad » 1°. Therefore 
one finds that, in order to aim to get an angular resolution of the order of A/3m ~ 1°, one needs 6 ~ 0.1° and 
5t ~ 2.5 fj,s. 

Note that, as long as St is smaller than the pixel transit time the error is dominated by the pixel size and 
not by timing. Moreover one might assume that, with a best fit, this results will scale as TV -1 / 2 . 



The angular resolution A/3|| as a function of 9, is shown in figure 27 Some numerical estimates are given in 
the section [Ol 



5.9 X M resolution 

The requirement on AAm is satisfied following the previous requirements. In fact the relation between slant 
depth X and the coordinate along the EAS, £, is, in the approximation ^ and neglecting the Earth curvature, 
is: 

' dX 

In 



T^X. (38) 



Therefore for 5 re f the A£ corresponds AXm ~ 35 g/cm 2 (I'm ~ 824g/cm 2 ) is larger than 0.5 km. Therefore, 
neglecting the geometrical and kinematical details of the EAS development, it is larger than half the pixel size 
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Figure 27: Angular resolution A/?y as a function of 9 in the reference conditions. Solid line: H = 700 km; dotted 
line: H — 400 km (the line for H — 1000 km is very close to the first one). We assumed #y ~ 1 A*s. 



projected at the Eartrj^] The required Xm resolution is therefore compatible with the pixel size. An improved 
resolution of Xm ~ 20 g/cm 2 might be reached by a suitable fitting procedure. 



5.10 Noise and background 

The RB in the wavelength range WR is estimated from the results of the experiments BABY [3"9~] . 
NIGHTGLOW gO] , Arizona- Airglow gT] and the Universitetsky-Tatiana microsatellite 02] . 

The noise level generated by all the parts of the experimental apparatus has to be kept well below the physical 
background level, due to the faintness of the air scintillation signal. 

The total RB rate intercepted by the apparatus (on the whole entrance pupil and full FoV) is: 



BAir sin 2 [7m] 



(39) 



A typical value for the total rate of RB hits detected by the apparatus is of the order of a few hundreds of 
MHz (see section 6.4 1. Therefore the intrinsic noise of the whole apparatus (all sources) is required to be less 



than a few GHz (over the whole Photo-Detector) in order to be negligible with respect to the RB. Note that this 
requirement should also include the stray-light coming from lack of light-tightness of the apparatus. 

In order to detect the faintest EAS reaching the required energy threshold it is mandatory to subtract online 
the RB as its rate is significant with respect to the air scintillation signal. Therefore a continuous monitoring of 
the average background on a pixel-by-pixel basis is, most likely, unavoidable to go low in energy in order to have 
a measure of the RB on the space-time scale comparable to the space-time scale of the EAS development. 

As a preliminary step for a space-based mission it is therefore necessary to carry on a detailed characterization 
of the RB, on the space-time scale of the EAS development, to improve our knowledge of it and improve 
background rejection [43] . 

Many experiments have measured this quantity. However in order to devise a method for online background 
subtraction it is mandatory to have a finer characterization of the space-time behavior of the RB on the space-time 
scale of the EAS development: Aa ~ 0.1° and At ~ /is. It is important to stress that the characterization should 
include measurements along different directions from nadir to cope with off- field RB. 



J For an orbital height H ~ 400 km and a pixel granularity of Aa ~ 0.1°, the pixel size projected on ground is ~ 0.8 km. 
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6 Some order of magnitude estimates 



6.1 Air scintillation signal 

Consider a typical hadron-induiced EAS with E = 10 19 eV, at 7 = 15°, with an ideal optics (eqH = cos W)- 
The signal time-integrated irradiance, dM / dA, reaching the apparatus is given in table [6] for different values 

of 9 and H (the other parameter values are the one listed in table [lj . 
Note the the irradiance scales as ~ (cosj/H) 2 . 



dM /dA 


ps [ph-m 2 ] 


9 = 30° 


9 = 50° 


9 = 70° 


H = 


400 km 


40 


70 


180 


H = 


700 km 


15 


25 


GO 


H = 


1000 km 


5 


10 


30 



Table 6: Time-integrated irradiance ——r of an EAS. 

dA 



The detection of such an EAS obviously requires a large entrance pupil of many squared meters (see next 
section). 

The typical angular length and the apparent time duration of the signal are given respectively in table [7] 
and GO 



£eas « [°] 


9 = 30° 


9 = 50° 


8 = 70° 


H = 400 km 


0.8 


1.7 


3.1 


H = 700 km 


0.4 


1.0 


1.8 


H = 1000 km 


0.3 


0.7 


1.2 



Table 7: Angular length £ E as of an EAS. 



Teas ~ [H 


9 = 30° 


9 = 50° 


9 = 70° 


H = 400 km 


60 


84 


102 


H = 700 km 


60 


84 


103 


H = 1000 km 


60 


84 


104 



Table 8: Apparent time duration Teas of an EAS. 



6.2 Requirements on the optical triggering efficacy 

Following the discussion in section [5~5| we can calculate a required overall photo-detection efficacy in the reference 
conditions 

£p D (7=15°)«^ph ( ^-J 

where we remember that 7V p h — 100. Using the values of the time-integrated irradiance given in the previous 
section we obtain £pd(7 = 15°) ~ 4 m 2 when observing at H — 700 km, £pd(7 = 15°) ~ 1.5 m 2 at H = 400 km 
and £ PD (7 = 15°) « 10 m 2 at H = 1000 km. 

If one conservatively assumes a safe overall (including filters, electronics, photo-sensor and filling factor) 



efficiency £prj ~ 0.1, we obtain an optical triggering efficacy, from the (111 




at H = 400 km 

4)( 7 = 15°) ps ^ w {40 m 2 ati? = 700 km . (40) 

at H = 1000 km 

The triggering optical efficacy at 7 = 0°, Ae' o (0°), i.e. the entrance pupil area, is not only the effective 
collection area but it is also an estimate of the physical area of the optics, actually an optimistic one, and 
consequently a requirement on the minimum size of the telescope. 
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Translated into minimum area and diameter for the optics we obtain, for an ideal optics in which £(7) = 
5*008 7, 

at // 100 km 

=£o(0°) « ^y.TJ ~ \ 42 m 2 at if = 700 km 

at i? = 1000 km 

(41) 

4.5 m at H = 400 km 




'As' (0°) 

Dnnu = -2\l — ^-^ ft < 7.5 m at H = 700 km 

12 m at H = 1000 km 

This result cannot be applied without accounting for the various sources of inefficiencies in a real optical 
system. Based on existing designs these inefficiencies can be assumed to reduce the amount of photons in the 
bucket depending on many factors. Therefore, the previous lower limit on D m [ n might rise by a factor of ^ 1, as 
confirmed by Monte-Carlo simulations [5J HJ [MJ |3S] . 



6.3 Granularity and angular resolution 

If we require a spatial granularity on Earth surface At < 1 km, we can fix the corresponding pixel granularity at 
different orbital heights: 

(0.10° atff = 400km 

A£ 

Att«-K 0.06° at H = 700 km . (42) 
[o.04° at #= 1000 km 

Using the values of table [7J one can easily see that at 9 = 50°, the EAS is about 17 pixels long; so there are 
about 6 photons per pixel. 



The pixel size on the FS is, from ( 18 1, S ~ Df/$= A£ / H sa 5 mm, using the minimum diameter given in (41 1 
and a //# ~ 0.5. The value of 5 doesn't depend on H since -Dmin oc H . 

With these numbers one can see that, in the reference conditions, the angular resolution perpendicular and 

parallel to the line of sight are both ~ 0.5°. The total angular resolution is then A/3 tot = J A/3j_ 2 + Af3\ 2 < 1°. 



6.4 The random background (RB) 

The total RB rate intercepted on the whole entrance pupil by the apparatus in the reference conditions, with 
7m = 20° (corresponding to 0.38 sr full FoV) is given in table [9] together with the corresponding number of 
pixels and the total RB rate detected per pixel (assuming a total photon efficiency, optics plus photo-detector, 
e T ~ 0.1). 

H [km] jV pix ggtot [THz] RB pi * [MHz] 
400 1.2-10 5 2^8 2T3 

700 3.5-10 5 7.6 2.2 

1000 7.8-10 5 19.0 2.4 



Table 9: Number of pixels, total RB rate (= B An sin 2 [7m]) and RB rate per pixel (= B An sin 2 [7 M ] e^/N). We 



used the values of A and Aa from (41) and (42 1 respectively 



This gives: one order of magnitude more RB than signal photons superimposed on the typical EAS (all 
space-time length) and roughly the same number of signal and RB photons near the EAS maximum. The 
acceptable background level also depends on the energy of the EAS. However to allow for background dependent 
observations implies a very precise knowledge of the Apparatus sensitivity as a function of the background level. 



7 Conclusions 

The design of a space-based apparatus for the detection of the EAS produced by UHECP is a challenging task, 
requiring a careful design and a preparation based on both preliminary studies and measurements and various 
technological demonstrators. 

Some sort of path-finder and/or technology demonstrator is certainly needed. 
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